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Abstract. Numerical studies of the interplanetary "shock overtaking mag- 
netic cloud (MC)" event are continued by a 2.5 dimensional magnetohydro- 
dynamic (MHD) model in heliospheric meridional plane. Interplanetary di- 
rect collision (DC) /oblique collision (OC) between an MC and a shock re- 
sults from their same/different initial propagation orientations. For radially 
erupted MC and shock in solar corona, the orientations are only determined 
respectively by their heliographic locations. OC is investigated in contrast 
with the results in DC [Xiong et al., 2006]. The shock front behaves as a smooth 
arc. The cannibalized part of MC is highly compressed by the shock front 
along its normal. As the shock propagates gradually into the preceding MC 
body, the most violent interaction is transferred sideways with an accompa- 
nying significant narrowing of the MC's angular width. The opposite deflec- 
tions of MC body and shock aphelion in OC occur simultaneously through 
the process of the shock penetrating the MC. After the shock's passage, the 
MC is restored to its oblate morphology. With the decrease of MC-shock com- 
mencement interval, the shock front at 1 AU traverses MC body and is re- 
sponsible for the same change trend of the latitude of the greatest geoeffec- 
tiveness of MC-shock compound. Regardless of shock orientation, shock pen- 
etration location regarding the maximum geoeffectiveness is right at MC core 
on the condition of very strong shock intensity. An appropriate angular dif- 
ference between the initial eruption of an MC and an overtaking shock leads 
to the maximum deflection of the MC body. The larger the shock intensity 
is, the greater is the deflection angle. The interaction of MCs with other dis- 
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turbances could be a cause of deflected propagation of interplanetary coro- 
nal mass ejection (ICME). 
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1. Introduction 

Interplanetary (IP) space is permeated by highly fluctuating solar wind with magnetic 
fleld frozen in its plasma [Parker, 1963]. The relatively quiet equilibrium of IP space is 
frequently interrupted by the solar disturbances, especially during solar maximum. Giant 
clouds of ionized gas with magnetic flux of 10^^ maxwell and plasma mass of 10^^ g, called 
coronal mass ejection (CME), are regularly emitted from the sun [Gosling, 1990; Webb et 
al., 1994]. IP CME (ICME) generally causes strong perturbation in the space environment 
as it passes by. Several models have already been applied in space weather forecasting, 
such as (1) HAF (Hakamada-Akasofu-Fry) [Fry et al, 2001, 2005]; (2) STOA (Shock 
Time of Arrival) [Smart and Shea, 1985]; (3) ISPM (Interplanetary Shock Propagation 
Model) [Smith and Dryer, 1990]; (4) an ensemble of HAF, STOA and ISPM models 
[Dryer et al., 2001, 2004]; (5) SWMF (Space Weather Modeling Framework) [Groth et al., 
2000; Gombosi et al., 2001; Toth et al., 2005]; (6) HHMS (Hybrid Heliospheric Modeling 
System) [Detman et al., 2006] , and so on. Great challenges are still faced to improve 
the prediction performance of space weather to satisfy the ever-increasing demands from 
human civilization [Baker, 2002]. 

Magnetic clouds (MCs) are an important subset of ICMEs, whose fraction decreases 
from ~ 100% (though with low statistics) at solar minimum to ~ 15% at solar maximum 
[Richardson and Cane, 2004, 2005]. Identifled by their characteristics including enhanced 
magnetic fleld, large and smooth rotation of magnetic fleld and low proton temperature 
[Burlaga et al., 1981], MCs have been the subject of increasingly intense study. The 
MCs with long interval of large southward magnetic fleld Bg are widely considered to 
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be the major IP origin of moderate to intense geomagnetic storms, especially during the 
solar maximum [Tsurutani, 1988; Gosling et al., 1991; Gonzalez et al., 1999] and, hence, 
play a crucial role in space weather prediction. An MC should probably be a curved 
loop-like structure with its feet connecting to the solar surface [Larson et al., 1997]. The 
force-free magnetic flux rope models have been proven to be very valuable to interpret 
in situ observations of MCs [Lundquist, 1950; Goldstein, 1983; Burlaga, 1988; Farrugia 
et al., 1993]. For the study of evolution of an individual MC during its anti-sunward 
propagation, many sophisticated models are developed based on these initial flux rope 
models: (1) Analytical models [Osherovich et al, 1993a, b, 1995; Hidalgo, 2003, 2005]; (2) 
Kinematic models [Riley and Crooker, 2004; Owens et al., 2006]; (3) Numerical models 
[Vandas et al., 1995, 1996, 1997, 2002; Groth et al., 2000; Odstrcil et al., 2002; Schmidt 
and Cargill, 2003; Manchester et al., 2004a, b]. Especially numerical simulations in (3) on 
a single MC have been exhaustive under the condition of various magnetic field strengths, 
axis orientations and speeds. 

ICME is not an absolutely self-isolated entity during IP propagation. It may interact 
with other solar transients (e.g., shock, ejecta) and heterogenous medium (e.g., corotat- 
ing interacting region). With less defined characteristics, some IP complex structures are 
reported recently, such as complex ejecta [Burlaga et al., 2002], multiple MCs [Wang et 
al., 2002a, 2003a], shock-penetrated MC [Wang et al., 2003b; Berdichevsky et al., 2005], 
non-pressure-balanced "MC boundary layer" associated with magnetic reconnection [Wei 
et al., 2003, 2006], ICME compressed by the following high-speed stream [Dal Lago et al., 
2006], and so on. Dynamical response and ensuing geoeffectiveness of these structures are 
directly associated with the interaction during their formation and evolution. Numerical 
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simulations have been applied to study most of the complex structures: e.g., the inter- 
action of a shock wave with an MC [Vandas et al., 1997; Odstrcil et al., 2003; Xiong et 
al., 2006], and the interaction of two MCs [Odstrcil et al., 2003; Gonzalez-Esparza et al., 
2004; Lugaz et al., 2005; Wang et al., 2005]. 

The observed "shock overtaking MC" events substantiate the likelihood of strong shock 
propagation in low (3 medium of MC plasma and, therefore, present a very interesting 
topic in IP dynamics. The evolution stages of MC-shock interaction within 1 AU are de- 
termined by MC and shock commencement interval in solar corona. They can be assorted 
into two categories: (1) shock still in MC (e.g. October 3-6 2000 and November 5-7 2001 
events [Wang et al., 2003b]); (2) shock ahead of MC after completely penetrating it (e.g. 
March 20-21 2003 event [Berdichevsky et al., 2005]). The idea that shock compression of 
the preexisting southward magnetic component can increase geoeffectiveness of the cor- 
responding Bg event has been proved in data analyses [Wang et al., 2003d]. Particularly, 
MC-shock compounds in category (1) cause highly intense geomagnetic storms [Wang et 
al., 2003b, c; Xiong et al., 2006]. Furthermore the geoeffectiveness variance of MC-shock 
compound with respect to the increasing depth of a shock entering a preceding MC was 
investigated in our previous study [Xiong et al., 2006, hereinafter referred to as paper 
1]. Both MC core and shock nose are radially erupted along heliospheric current sheet 
(HCS) in paper 1; however, the above-mentioned specific MC-shock events [Wang et al., 
2003b; Berdichevsky et al., 2005] were all identified such that the shock flank sweeps the 
preceding MC body. IP direct collision (DC)/oblique collision (OC) of an MC and a shock 
results from their same/different initial propagation orientation. For radially erupted MC 
and shock in solar corona, the orientations are only determined respectively by the heli- 
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ographic locations of MC core and shock nose. Because the probabihty of MC core and 
shock nose radially launching from the same heliographic location is very rare and shock 
front extends over a wide angular span in IP medium, it is meaningful to study the role of 
shock orientation relative to a preceding MC propagation. DC in paper 1 is here modified 
to be OC for MC-shock interaction. The shock in DC/OC is correspondingly named as 
"central" / "non-central" shock. Moreover DC/OC is likely to be the IP interaction of two 
radially propagating disturbances from the same/different solar activity regions. 

Section 2 presents a brief description of numerical magnetohydrodynamic (MHD) model. 
Section 3 discusses the dynamical evolution of MC-shock OC. Section 4 analyzes the 
ensuing geoeffectiveness of MC-shock compound. Section 5 describes the dependence of 
shock-induced MC deflection on shock orientation and intensity. Section 6 summaries the 
conclusions. 

2. Numerical MHD Model 

The detailed description of the numerical model, including numerical scheme, compu- 
tational mesh layout, prescription of the ambient solar wind and preceding MC, is given 
in paper 1. Only the shock introduction among input parameters of numerical model is 
modified to simulate OC of MC-shock interaction in contrast with DC in paper 1. 

An incidental fast shock, which is radially launched from the inner boundary, is pre- 
scribed by several parameters: its emergence time tso, the latitude of its nose 9sc, the 
latitudinal width of its flank A^^, the maximum shock speed within its front Vs, the du- 
ration of growth, maintenance and recovery phases {tgi, ts2, tgs). Some parameters are 
fixed in all simulation cases of paper 1 and here, i.e. A^^ = 6°,tsi = 0.3 hour, = 
1 hour,t53 = 0.3 hour. The remaining parameters (tsO) ^sc? Vg) are independently chosen 
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to mimic different conditions of IP MC-shock interaction, tso is used to separate the MC 
and shock initiahzation in time for reproducing different evolutionary stages of MC-shock 
compound at 1 AU. 9sc designates emergence orientation of shock nose relative to pre- 
vious MC propagation. Since the preceding MC emerges from the heliospheric equator, 
Bsc = 0° and Ogc 7^ 0°, corresponding to the introduction of "central" and "non-central" 
shock, determine MC-shock DC and OC in IP space respectively. Vg describes the inten- 
sity of MC-shock interaction to some extent. All introduced shocks in our simulation are 
strong enough to be faster than the local magnetosonic speed at all time and, therefore, 
to prevent weak shock dissipation in MC medium. 

3. Dynamics of MC-shock Interaction 

All fifty simulation cases are assorted into five groups in Table 1. Groups of individual 
MC (IM), direct collision (DC), oblique collision (OC), shock orientation dependence 
(SOD), and shock intensity dependence (SID) are studied respectively, where Groups IM 
and DC have been addressed in detail in paper 1. Case Pi is shared by Groups DC 
and SOD, and Case P2 shared by Groups OC, SOD and SID. With the identical Vg of 
1630 kms~^ and variable tso from 3 hours to 41 hour. Groups DC and OC only differ in 
9sc for comparative study. By modifying 9sc from 0° to 10°, "central" shock in DC is 
directed to be "non-central" one in OC. Further, the parametric studies of 9sc from 0° to 
45° in Group SOD and Vs from 947 kms~^ to 3173 kms~^ in Group SID are explored as a 
supplement to Groups DC and OC. Cases Bi and B2 with tgo = 41 hours, and Cases Ci 
and C2 with t^o = 10 hours are typical examples of MC-shock interaction in categories 1 
and 2 referred in Section 1. 
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3.1. Case B2 

The process of MC-shock interaction of Case B2 is visualized in Figure 1. Under each 
image are two corresponding radial profiles by cutting right through 0° (noted by Lat. = 
0°) and southern 4.5° (white dashed lines in the images, noted by Lat. = 4.5°S) away from 
the equator. The magnitude of magnetic field in radial profile is given by subtracting its 
corresponding initial value of ambient equilibrium. The body of MC is identified to be 
enclosed by a white solid line in the images and between two dotted lines in attached 
profiles. Magnetic field configuration is superimposed upon the images. The incidental 
shock aphelion arrives at 90Rs (along Lat. = 4.5°S) in 50.4 hours meanwhile the MC 
core arrives at IGORg (along Lat. = 0°), shown in Figure 1(a), (d) and (g). Impending 
collision can be pregnant from large radial speed difference between the preceding MC and 
the following shock, as indicated by radial bulk flow speed Vr of 830 kms~^ at shock front 
and 540 kms~^ at MC head from the profile of Lat. = 4.5°S (Figure 1(d)). Though the 
latitudinal span of its flank is 6° initially at inner boundary, the shock extends up to 40° 
quickly due to its very strong intensity, until it emerges into IP medium completely. The 
traverse of shock front across the equator leads to significant HCS warping seen clearly 
in Figure 1(b), which is consistent with previous results [Smith et al., 1998; Hu and Jia, 
2001]. As shock emergence orientation is redirected, the morphology of IP shock changes 
from a concave (Figure 3(e) in paper 1) to a smooth arc (Figure 1(e) here). As a result, 
MC-shock interaction consequently changes from DC to OC. The shock just catches up 
with the inner boundary of MC at 66.9 hours (Figure 1(b), (e) and (h)). Due to strong 
magnetic field and low (3 plasma, the radial characteristic speed of fast mode wave Cf 
of the MC is abnormally high at 1 AU with 200 kms~^ in maximum at MC core and 
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100 kms~^ in minimum at MC boundary. The rare chance of shock survival in an MC 
medium explains why only a few "shock overtaking MC" events are observed in IP space. 
Across the tangent point between inner MC boundary and shock front exists a quite 
sharp slope of Vr, as clearly seen along Lat. = 4.5°S. MC-shock interaction begins from 
this tangent point at 66.9 hours. Once a slow MC is within the very large latitudinal 
span of the overtaking shock front, it will be swept by the shock and, from then on, the 
evolution of MC and shock will be coupled with each other. The overwhelming shock 
significantly distorts MC morphology at 85.4 hours (Figure 1(c), (f) and (i)). Namely, 
the originally curved magnetic field lines become very fiat. The collision is more violent 
along Lat. = 4.5°S. A sharp discontinuity is conspicuously formed in the rear part of MC 
with B — B\t=o = 25 nT, Vr = 620 kms~^, and cj = 260 kms~^ in maximum within highly 
compressed region. 

3.2. Case Ca 

In Case C2, an earlier shock emergence (t^o = 10 hours) allows the incidental shock to 
ultimately penetrate the MC body within the solar-terrestrial heliospheric range. Only 
the evolution of Vr is given in Figure 2 to show the concerned MC-shock complex structure. 
Though an MC generally behaves like a rigid body with a little elasticity, magnetic field 
lines of the simulated MC appear to be too vulnerable to be easily deformed in the face 
of an overwhelming shock. The shock is radially emitted with the strongest intensity at 
front nose. Hence shock front behaves as an oblique curve relative to heliospheric equator 
due to the propagation speed difference from shock nose to edge flank. The MC is highly 
compressed by the shock along its normal. The shock front looks like a smooth arc in 
MC medium. As it propagates gradually into the preceding MC body, the most violent 
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interaction is transferred sideways (heliolatitudinally in the present study). Due to net 
shock-input angular momentum during MC-shock OC, the MC core starts to deflect away 
from initial shock orientation when the shock enters MC core, as seen in the contrast of 
Figure 2(b) and (c). The overall MC body is also deflected to the north. The global MC 
body deflection is quantified by the deflection angle of its core. Once the shock completely 
penetrates the MC, the grip of shock force on the MC is substantially relaxed, and the 
MC is restored to the roughly ellipse morphology by its field line elasticity. Meanwhile, 
the MC loses its angular speed component by the relative difference between the radial 
ambient flow and the speed's value at the MC boundary and, then propagates radially 
along the deflected angle. The incidental shock is also simultaneously deviated with its 
aphelion in the opposite direction, until it finally merges with the MC-driven shock into 
a compound one. The bend of interplanetary magnetic field (IMF) lines is obvious near 
the south of MC boundary, seen from Figure 2(c). 

Figure 3 shows the comparison among Cases A, Ci, and C2 about time- dependent 
parameters: (a) radial distance of MC core r^, (b) MC radial span Sr, (c) MC angular 
span S9, (d) MC cross section area A, and (e) MC core deflection angle D9m, where the 
solid, dashed and dashed-dotted curves denote Cases A, C2 and Ci respectively, and the 
three vertical delimiting lines (dotted, dashed and dotted) from left to right correspond 
to the occasion of shock encountering MC tail, core and head respectively. The MC in 
Case C2 is largely compressed by the shock, beginning from 13 hours. The dependence 
of the compression of MC geometry on shock orientation is illustrated by the comparison 
in Figure 3(b)-(d). Sr is larger while S9 is smaller for Case C2 in Group OC. Though 
S9 is little affected in Case Ci when shock front is in MC body (13 hrs < t < 33 hrs). 
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it is significantly narrowed in Case C2. And the MC cross section area A in Case C2, 
which represents the overall influence of shock compression due to integration of factor 
Sr and S9, is a bit larger than that in Case Ci. Starting from being encountered by 
the following shock, MC core deflects up to —4.5° until shock front reaches MC head, as 
seen in Figure 3(e). Though total deflection angle of MC core (—4.5°) amounts to three 
computational grids of latitudinal spacing 1.5°, MC deflection, we think, is indeed physical 
solution. Due to rough subcell resolution in numerical computation, MC core deflection 
behaves as a false discrete quantum-like transition instead of a realistic smooth one. But 
it does not distort the fundamental physical characteristics in numerical simulation. 

3.3. Multi-Cases Comparison 

The propagation of MC-shock structure toward the earth can be detected by Ll-orbiting 
spacecraft, which perform the sentinel duty in space weather alarm system. The montage 
of the evolution of MC-shock compound at LI under three typical circumstances is visu- 
alized in Figure 4, where (a)-(c) correspond to Case Ri from Group DC and Cases Q2 
and R2 from Group OC. Though the farthest radial distances of shock front in the north 
and south of the equator are almost identical in Cases Ri and Q2, the shock intensity in 
the south in Case Q2 is apparently stronger than its north counterpart. With a smaller 
emergence interval, the shock in Case R2 merges completely with the MC-driven shock 
into a compound one and moves faster in the south by contrast of Figure 4(b) and (c). 
Moreover the asymmetry of compound shock front with respect to heliospheric equator 
occurs when the shock erupts sideways relative to the MC propagation. The final MC 
propagation is slightly deviated from heliospheric equator to northern 4.5° after being 
ultimately penetrated by the shock, as seen from Figure 4(b) and (c). The succedent 
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high speed flow right after the inner boundary of preceding MC in Group DC, mentioned 
in paper 1, does not exist in corresponding Group OC, which can be seen from contrast 
between Figure 4(a), (b), and (c). The shock front with 9^^ 7^ 0° has the obhque nor- 
mal relative to the preceding MC propagation, so the disturbance of speed enhancement 
downstream of shock front in Group OC can completely bypass or penetrate the obstacle 
of MC body and merge with the MC-driven shock. 

4. Geoeffectiveness Studies 

The southward magnetic flux within the MC is located in its rear part. The geomagnetic 
effect of simulated Bg event is quantified by Dst index. The in-situ measurements by 
a hypothetic spacecraft at LI are inputted to Burton formula [Burton et al., 1975] to 
calculate Dst, as applied by Wang et al. [2003c]; Xiong et al. [2006]. 

Near-HCS latitudinal dependence of Dst index in Cases A, B2, and C2 is plotted in 
Figure 5. The positive and negative latitudes are referred to southern and northern semi- 
heliosphere. With the MC core marked by A and MC boundary by o, the solid, dashed, 
and dashed-dotted lines denote Cases A, B2, and C2 respectively. Geomagnetic storm 
has been obviously aggravated by shock overtaking MC. The minimum Dst are found 
to be -103 nT, -168 nT, and -140 nT in Cases A, B2, and C2, respectively. Cases B2 
and C2 are discussed one by one against Case A. Firstly, geomagnetic storm in Case A 
is largely enhanced in Case B2 within the latitudinal span influenced by the shock. The 
minimum Dst occurs at 3° rather than 0° (the latitude of MC core passage), because 
the former undergoes more violent compression. The geoeffectiveness remain unchanged 
within Lat. < —5°. The asymmetry of shock propagation with respect to heliospheric 
equator leads to subsequent asymmetry of geoeffectiveness of the MC-shock compound. 
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Secondly, in Case C2 the concave of the latitudinal distribution of Dst is shifted 4.5° to 
the north. The MC deflection is caused by "non-central" shock penetrating MC body, as 
interpreted in section 3.2. As a result, the southward passing magnetic flux decreases due 
to the northward deflection of MC, and the IMF bend south of the equator due to shock 
passage, seen from Figure 2(c), which are responsible for the increased and decreased 
Dst in 2.3° < Lat. < 9.4° and 9.4° < Lat. < 15° respectively, comparing with Case A. 
Therefore, as shock front propagates from the south (Case B2) to the north (Case C2) in 
MC medium, the latitude of minimum Dst consequently moves in the same direction. 

All MC-shock interaction cases of Group OC are integrated to study further the de- 
pendence of Dst index on the penetration depth dost of shock overtaking MC. d£,st is 
defined as the radial distance between shock front and MC inner boundary along sun-MC 
core. Three in-situ observations in time sequence at LI along heliospheric equator and 
±4.5° aside are synthetically analyzed in Figure 6, where the three vertical delimiting 
lines (dotted, dashed and dotted) from left to right correspond to the cases of shock en- 
countering the tail, the core and the front of MC at LI, respectively. From top to bottom 
are plotted (a) MC-shock emergence interval, noted by Dt, (b) the Dst index, (c) the 
minimum of dawn-dusk electric field VB^., noted by Mm.iV Bz)., (d) the interval between 
the commencement of VB^ < —0.5 mV/m and the corresponding Dst minimum, noted by 
At, (e) the minimum of southward magnetic component Bg, noted by Min. (Bs), and (f) 
the maximum of magnetic field magnitude Max.(i?), respectively. The solid, dashed and 
dashed-dotted lines in Figure 6(b)-(f) correspond to the observations at Lat.=0°, 4.5°S 
and 4.5°N, respectively. The separate MC and shock events are coupled together when 
Dt < 50 hours. The shock penetrates into the preceding MC more deeply with shorter 
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Dt. Min. (Bs) and Mm.{V Bz) decline dramatically along Lat. = 4.5°S as dost increases 
from to lORs, because the first tangent point between MC boundary and shock front 
is very near 4.5°S. Dst decreases monotonically within ORg < d^gt < 23. 5/?^ until shock 
front reaches MC core. Once the shock front exceeds the MC core {dost > 23.5i?s), the 
latter begins to deflect northward. Moreover when dost > 23.5/2^, the greatest compres- 
sion region by the shock front is within the MC anterior part or the MC-driven sheath, 
where magnetic field is northward and, hence, contributes little to geoeffectiveness. So 
the mitigated geoeffectiveness along 0°, 4.5°S and aggravated geoeffectiveness along 4.5°N 
coexist, as seen from 23. 5/?^ < dost < 44. 5/?^ in Figure 6(b). 

Based on the analyses of Figures 5 and 6, MC deflection by MC-shock OC plays a crucial 
role in geomagnetic storms. The minimum Dst and its corresponding latitude among Dst 
latitudinal distribution for every case of Group OC are assembled in Figure 7. With a 
given Dt, there exists a latitude where geoeffectiveness reaches its maximum (Figure 7 
(a)). This specific Dst value is plotted as dashed line in Figure 7(b). The latitudinal 
distribution of individual MC event (Case A), serving as a background in contrast, is also 
plotted as solid line in Figure 7(b). The relative ratio of geoeffectiveness enhancement 
by the shock is presented in Figure 7(c) to quantify two curves difference in Figure 7(b). 
As Dt decreases from 48 hours to 3 hours, the latitude of maximum geoeffectiveness 
firstly remains constant with decreased Dst from -115 nT to -180 nT, enhanced ratio 
from 20% to 91%, then monotonically changes from 3° to —4.5° with gradually subdued 
geoeffectiveness, finally remains constant again with further increased Dst from -130 nT 
to -115 nT, decreased ratio from 50% to 30%. The minimum Dst (-185 nT) occurs at 
2.3° when the shock front enters MC core right at 1 AU. In contrast with paper 1, the 
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maximum geoeffectiveness of MC-shock interaction in Group DC is the same as that in 
Group OC despite occurrence at different heholatitudes. 

5. MC and Shock Deflections 

IP MC deflection mentioned in Section 3.2 is a key parameter for solar-terrestrial trans- 
portation process, because it concerns the preexisting condition of geomagnetic storms 
- whether an MC could encounter the earth. In order to explore reliance of MC core 
deflection angle on shock orientation and intensity, the results of Groups SOD and SID 
are illustrated in Figure 8. Because MC core continuously deflects on the condition of 
shock front being in MC medium, seen from Figure 3(e), all tgo in Groups SOD and SID 
are chosen to be 10 hours to have MC completely penetrated for obtaining final invariant 
angular displacement of MC core 06^- Dst in Figure 8 refers to the geoeffectiveness at 
certain latitude of passage of deflected MC core. Firstly for Group SOD with different 
6sc, two factors affect D9m- (1) Osc 7^ is a premise of MC core deflection; D9m = 
corresponds to = 0. (2) As increases, shock flank section encountered by MC body 
is further away from shock nose and, hence weaker. The absolute value of deflection angle 
tends to be smaller due to the weakening of MC-shock collision. The maximum deflection 
of MC core {DOm = -4.5°) occurs at certain Osc (10° < Osc < 15°). Meanwhile Dst in- 
creases monotonically as a function of Ogc, up to the value of corresponding individual MC 
event. Secondly for Group SID with different Vg, both D6m and Dst decrease steadily as 
Vs increases. Moreover, the slopes of two curves in Figure 8(c) and (d) decrease steadily, 
very abrupt in Vg = 1000 km/s and become nearly horizontal when Vg > 3000 km/s. This 
saturation effect on D9m and Dst is caused by the concurring deflection of shock aphelion 
opposite to that of MC core mentioned in Section 3.2. So the divergent trend of deflection 



DRAFT 



April 5, 2009, 12:26pm 



DRAFT 



XIONG ET AL.: MC-SHOCK INTERACTION AND ITS GEOEFFECTIVENESS 2 X - 17 

angle between the MC body and the shock aphehon counteracts, more or less, the effect 
of increasing shock speed Vg on MC-shock collision. 

The finding of MC deflection due to interaction with a shock is further discussed through 
comparison with other relevant models. (1) Vandas et al. [1996] proposed that an MC 
deflects during the propagation through IP medium with unipolar IMF. Magnetic recon- 
nection between IMF and inherent MC field across one side of MC boundary causes the 
angular force unbalance and, hence, leads to angular deflection. The MC continuously 
deflects through IP space. The role of magnetic helicity is responsible for deflection mech- 
anism [Vandas et al., 1996]. However, such deflection needs to be verified further, as 
the reconnection should not be so significant in the IP medium with low f3] (2) Wang et 
al. [2004] suggested that CMEs could be deflected as largely as several tens degrees in 
the propagation under the effects of background solar wind and spiral IMF. CME deflects 
from its onset until accelerated or decelerated to background solar wind, which is expected 
to be done within several tens solar radii [Wang et al., 2006b]. It can well interpret the 
observation fact of east-west asymmetry of solar source distribution of earth-encountered 
halo CMEs [Wang et al., 2002b] and why some eastern limb CMEs encountered the earth 
[Zhang et al., 2003] and some disk CMEs missed the earth [e.g., Schwenn et al., 2005; 
Wang et al., 2006a]; (3) Our model here gives that MC deflection only happens during 
the process of shock front penetrating MC body. The effect of shock pushing MC aside 
leads to the deviation of MC by several degrees at the most; (4) We conjecture that 
interaction between ICMEs may also be a cause of ICME deflection, and the deflection 
angle could be up to tens degrees, larger than that in (3). The propagation trajectory of 
CMEs mentioned above is deflected from an initial straight line in the IP medium. Both 
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deflections in (1) and (2) are caused by interaction between ambient solar wind and IP 
disturbance. In contrary, the deflection in (3) and (4) are ascribed to interaction between 
different IP disturbances, i.e. the collision between MC-shock or MC-MC. It may expect 
a significant effect on the possibility of CME hitting the earth in (1), (2), and (4), whereas 
the effect in (3) may be negligible because of the small deflection angle. 

The deflection of shock aphelion in IP medium is a key factor in the near-earth prediction 
of shock arrival time. Hu [1998]; Hu and Jia [2001] stated that the deflection of shock 
aphelion results from joint effects of spiral IMF and heterogenous medium consisting of 
fast and slow solar wind. The deflection is also found here in OC of MC-shock. Starting 
from shock passage through MC medium, shock aphelion deflects toward the contrary 
trend of MC deflection until the shock totally merges with the MC-driven shock. The 
flnal shock aphelion as well as front morphology are distinct from those of isolated shock 
event. Both MC and shock undergo signiflcant modification during the process of their 
collision. 

6. Concluding Remarks and Discussions 

For further understanding of the IP "shock overtaking MC" events [Wang et al., 2003b; 
Berdichevsky et al., 2005], the investigation of MC-shock interaction and consequent geo- 
effectiveness in paper 1 is continued by a 2.5-dimensional ideal MHD numerical model. 
The simulations find that shock eruption orientation relative to preceding MC propagation 
plays a crucial role in MC-shock interaction. 

Firstly, MC-shock dynamical interaction is modeled. In order to reveal the effect of the 
shock orientation relative to preceding MC propagation, DC in paper 1 is here modified to 
be OC for MC-shock interaction under the condition of the same shock speed. The results 
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show that the shock front in MC-shock OC behaves as a smooth arc in MC medium. The 
cannibahzed part of MC is highly compressed by the shock along its normal. As the shock 
propagates gradually into the preceding MC body, the most violent interaction is trans- 
ferred sideways (in terms of heliolatitude) with an accompanying significant narrowing 
of the MC's angular width. The opposite deflections of MC body and incidental shock 
aphelion concur during the process of shock penetrating MC. MC deflection ends when 
the shock approaches MC head; Shock deflection stops when the shock completely merges 
with MC-driven shock. After shock passage the MC is restored to oblate morphology. 
The high speed flow right after MC inner boundary mentioned in paper 1 does not exist 
here on the condition of non-uniform orientation of initial MC and shock eruption. 

Secondly, the geoeffectiveness of MC-shock OC is studied. Geoeffectiveness of an indi- 
vidual MC is largely enhanced by an incidental "non-central" shock. With the decrease 
of MC-shock commencement interval, shock front at 1 AU traverses MC body and is 
responsible for the same change trend of the latitude of the greatest geoeffectiveness of 
MC-shock compound. Among all cases with penetrating shock at various stages, the 
maximum geoeffectiveness occurs when the shock enters MC core right at 1 AU. Wang 
et al. [2003c] suggested that the maximum geomagnetic storm be caused by shock pene- 
trating MC at a certain depth, and the stronger the incident shock is, the deeper is the 
position. Based on our numerical model, Wang's conclusion of shock penetration depth 
regarding the maximum geoeffectiveness [Wang et al., 2003c] may be supplemented that 
shock position is right at MC core on the condition of very strong shock. 

Thirdly, the reliance of MC deflection on shock orientation and intensity is explored. 
The angular displacements of MC body and shock aphelion are ascribed to MC-shock OC. 
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An appropriate angular difference between the initial eruption of an MC and an overtaking 
shock leads to the maximum deflection of the MC body. The larger the shock intensity is, 
the greater is the deflection angle. The interaction of MCs with other disturbances could 
be a cause of ICME's deflected propagation. 
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Figure Captions 

Figure 1 The evolution of shock overtaking MC for Case B2, with (a)-(c) magnetic field 
magnitude (d)-(f) radial flow speed Vr-, and (g)-(i) radial characteristic speed of fast 
mode Cf. Below each image are two additional radial profiles along Lat.= 0° and 4.5°S. 
Note: radial profile of B is plotted by subtracting initial ambient value -B|f=o- The white 
solid line in each image denotes the MC boundary. Solid and dashed lines at each profile 
denote MC core and boundary. 

Figure 2 The evolution of shock overtaking MC for Case C2 with radial flow speed v^-. 
Only part of domain is adaptively plotted to highlight MC. 

Figure 3 The time dependence of MC parameters: (a) radial distance of MC core r^, 
(b) MC radial span Sr , (c) MC angular span SO, (d) MC cross section area A, and 
(e) MC core deflection angle Ddm- The solid, dashed and dashed-dotted curves denote 
individual MC event (Case A), MC-shock events (Case C2, Ci). Three vertical delimiting 
lines (dotted, dashed and dotted) from left to right correspond to the occasion of shock 
encountering MC tail, core and head respectively. 

Figure 4 The montage of radial flow speed Vj. for the evolution of MC-shock compound 
at LI under three conditions. 

Figure 5 The comparison of latitudinal distribution of Dst index among individual MC 
event (Case A) and MC-shock events (Cases B2 and C2). The solid, dashed, and dashed- 
dotted lines denote Case A, B2, C2 respectively, with the mark A, o for the passage of 
MC core and boundary. The positive and negative latitude are referred to southern and 
northern semi-heliosphere. 
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Figure 6 The parameter variances of MC-related geoeffectiveness as a function of dnst in 
Group OC. Here dost refers to radial distance between shock front and MC inner boundary 
along sun-MC core. From left to right, three vertical lines (dotted, dashed, dotted) denote 
the occasions of shock just reaching MC tail, core, and front at LI respectively. The mark 
o and A denote corresponding results of Case B2 and C2. (a) Dt, MC-shock emergence 
interval, (b) Dst index, (c) Min.(yi?^), the minimum of dawn-dusk electric field VB^, (d) 
At, the interval between the commencement of VBz < —0.5 mV/m and the corresponding 
Dst minimum, (e) Min.(Ss), the minimum of southward magnetic component, and (f) 
Max. (5), the maximum of magnetic magnitude. Solid, dashed and dashed-dotted lines 
in (b) to (f) correspond to observations along Lat.=0°, 4.5°S and 4.5°N respectively. 

Figure 1 The response of the latitude of maximum geoeffectiveness and accompanying 
Dst (dashed line in (b)) as the change of MC-shock interval Dt in Group OC. The lat- 
itudinal distribution of individual MC event (Case A) is denoted in solid line of (b) as 
background. The relative ratio of geoeffectiveness enhancement (c) by the shock is derived 
from two curves difference of (b). 

Figure 8 The dependence of MC core deflection angle DOm and Dst at the specific 
latitude accompanying MC core passage, on shock eruption orientation (Group SOD) 
and speed Vg (Group SID). The horizontal dashed lines in (b) and (d) denote corresponding 
Dst of individual MC event (Case A). 
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Table 1. Assortment of simulation cases of individual MC and MC-shock interaction 
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Figure 1. 
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Figure 2. 
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Figure 3. 
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Figure 4. 
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Figure 5. 
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Figure 6. 
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Figure 8. 



DRAFT 



April 5, 2009, 12:26pm 



DRAFT 



